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The light variability of the helium strong star HD 37776 as a result 
of its inhomogeneous elemental surface distribution 
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Abstract. We simulate light curves of the helium strong chemically peculiar star HD 37776 assuming that the observed periodic 
light variations originate as a result of inhomogeneous horizontal distribution of chemical elements on the surface of a rotating 
star. We show that chemical peculiarity influences the monochromatic radiative flux, mainly due to bound-free processes. Using 
the model of the distribution of silicon and helium on HD 37776 surface, derived from spectroscopy, we calculate a photometric 
map of the surface and consequently the uvhy light curves of this star Basically, the predicted light curves agree in shape and 
amplitude with the observed ones. We conclude that the basic properties of variability of this helium strong chemically peculiar 
star can be understood in terms of the model of spots with peculiar chemical composition. 
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1. Introduction 

Periodic light variations are a common feature observed in the 
magnetic chemically peculiar (CP) stars. They can be inter- 
preted by the presence of "photometric" spots causing an un- 
even distribution of the emergent flux on the surface of rotating 
stars. Likely, these hypothetical photometric spots are some- 
how related to the "spectroscopic" spots of a peculiar abun- 
dance of various chemical elements. Within this paradigm, the 
departures of the energy distribution in spectra of CP stars from 
the energy distribution of normal stars of the same spectral type 
are as a rule explained as a consequence of line blanketing by 
plenty of lines originating in the spectroscopic spots with pe- 
culiar chemical co mposi t ion and the flux redistribut i on ind uced 
by diem (Peterson! ! 1975 lMolnaij[l973l) . lLanz et all d 19961) and 
Krticka et al.i (i2004 suggested that bound-free transitions may 



play a more important role in the blanketing, redistribution of 
the emergent flux, and the light variability. 

However, more effects may influence the spectral energy 
distribution and consequently may be invoked to explain the 
observed light variability. These effects may be connected with 
the influence of complex surface magnetic fields. IStaudd d 1 9721) 
and Trasco (1972) showed that the influence of the magnetic 
pressure on the stellar atmosphere structure migh t cause ob- 
servable light variations (see also Stepieiil 1978 : Carpenter 
Ll985; LeBlanc et al. 1994). The effect of the magnetic field 
on the atmosphere hydrostati c equilibrium was recently stud- 
ied by IValvavin et al.l (120041) . The influence of the Zeeman 
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effect on th e line blanketing may modify the emergent ra - 
diative flux dKochukhov et al.l 120051: iKhan & Shulvakll2006l) . 
Consequently, there are many different mechanisms that might 
result in the appearance of photometric spots on the stellar sur- 
face. 

Contrary to these photometric spots, the location and size 
(e.g., maps of distribution of individual chemical elements on 
the stellar surface) of the spectroscopic spots can be determined 
by Doppler imaging techniques (e.g., Khokhlovaet al. 1 9971) . 
However, it is not clear whether and how the photometric spots 
are related to the spectroscopic ones. 

Another explanation of the light varia bility of hot stars 
wit h strong surface magne tic fields exists. iNakaiimal d 1985b 
and ' Smith & Grootd d200l') proposed that the observed light 
variations may be caused by light absorption in rotating cir- 
cu mstellar clouds. Recent ly, this model was successfully used 
bv iTownsend et al. I d2005 l) for prediction of the light curve of 
the He-rich chemically peculiar star cr Ori E. The properties 
of these circumstellar clouds were inferred using the rigi dly 
rotating magnetosphere model dTownsend & Owo cki' 2005), in 
which the material blown from the star by the radiatively driven 
wind accumulates i n the potential minim a along individual 
magnetic field lines. iTownsend et al.l ( 2005 ) used free parame- 
ters characterizing an amount of absorption in the clouds to pre- 
dict the light curve of cr Ori E, as it was not known for how long 
the accumulation had proceeded. Consequently, they were able 
to reproduce the shape of the light curve, but not the amplitude 
itself. Moreover, cr Ori E shows hydrogen emission lines that 
indicate the presence of circumstellar environment. Last but not 
least, the light curve of cr Ori E is exceptional when compared 
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with the Hght curves of other CP stars. However, many CP stars 
showing Hght variabihty are hkely too cool to have sufficiently 
strong winds enabling formation of circumstellar clouds. Thus, 
the light variability of other CP stars may be caused by other 
mechanisms. 

The fact that CP stars in many cases show not only pecu- 
liar chemical composition, but also a variable spectrum - in- 
dicating the uneven distribution of chemical elements on their 
surfaces - tempts one to attribute the light variability to this 
uneven surface distribution of elements. This uneven distribu- 
tion, through the line blanketing and/or the bound-free absorp- 
tion, determines the emergent flux depending on the location 
on the stellar surface and may cause the variability of a rotating 
star. However, to our knowledge, a light curve of a CP star had 
never been calculated ab initio from the known distribution of 
chemical elements as derived by Doppler i maging techniques. 



The in itial step in this direction was done bv lKrivosheina et al 



imd). who were able to successfully reproduce the observed 
vby light curve of CU Vir modifying the temperature gradient 
in the regions of silicon overabundance. 

For our simulations of the light curves of chemically pe- 
culiar stars we selected HD 37776. This star r esiding in 
the reflection/ern ission nebula IC 431 ( van den Be rgh 196^ 
Finkbeinen 120031) belongs to well-studied helium-strong CP 
stars. It has a very strong complex sur face magnetic field 
with a significant quadrupolar compo nent jBorra & Landstreet 
[Thompson & Lands^ eet The observed periodic 

variat ions of helium lines ( Nissea, 19761 IPederse n & Thomse: 
19771) and lines of some lighter elements (fshore & Brow 
1990 ) are interpreted as a consequence of rotation of a star with 
uneven elemental surface distribution. Based on this model, 
surface distributions and abun dances of various elemen t s were 
derived from s p ectroscopy jBo hlender & Landstree3 1988 : 
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Veto et al.l I199U iKhokhlova et al. ,2000,). Consequently, fliis 



star seemed to us to be an appealing target for modelling of 
the light variations. 

In this work we present computed light curves of the He- 
strong CP star HD 37776. We used the surface distributions and 



abundances of various elements derived by Khokhlova et al 
(I2OOQ) . As silicon and helium create areas with large overabun- 
dances, and moreover, as it is known that silicon may influ- 
ence the observed spectral energy distribution of silicon-rich 



19961: iKrtickaet al 



CP st ars (lArtru &Lanzlll987l: iLanz et al. 
20041) . we concentrated on these two elements in our calcu- 



lations. The role of silicon may be more general because this 
element is overabundant in vast m ajority of light variable CP 
stars (including cool CP stars, e.g., Adelmanlll973 ). 



2. Basic assumptions 

2.1. Atmosphere parameters 

Table [T] lists the atmosphere parameters of HD 37776 
ado pted in this study. We adopted parameters derived 

|g 



by iGroote & KaufmannI (119811) which a re appropriate for 



HD 37776 spectral type (lHarmaneclll988 ) Even though the 
value \ogg - 4.5 used by Khokhlova et al. was higher, 

the value of the surface gravity adopted does not have a signifi- 



cant effect on the prediction of the light variability (see Sect.|5]l. 
The phases for light curves were computed according to the 
linear e phemeris derived by AdelmanI (11997 "). which was used 
also by Khokhlova et al. ( 200ol) for surface mapping. 



Table 1. Steflar parameters of HD 37776 



T^fi 22 000 K 

logg(CGS) 4.0 
Inclination ; 45° 
Abundance -5.4 < [He/H] < 2.3 
-2.4 < [Si/H] < 2.6 
Viurb 2kms"' 



Groote & K aufmann (1981) 
Groote & Kaufmann (1981) 
Khokhlova et al. (2000) 
Khokhlova et al. (2000) 



For our stud y we adopted the surface distribution of He and 
Si as derived by IKhokhlova et al. I (12000 *. see Fig. [B- Here the 
abundance is defined relatively to the hydrogen, i.e., [A/H] - 
logiNA/Nu) - log(A^A/Mi)0, where the arguments of the loga- 
rithms are number densities of elements A and H on the stel- 
lar surface, and on the solar surface respecti vely. Taking into 



accoun t that oxygen and iron as derived by IKhokhlova et al 



(I2OOOI) are largely underabundant on the entire surface, we did 
not consider these elements as sources of the light variability. 

2.2. Model atmospheres and photometric colours 

We use the standa r d co d e TLUSTY dHubenvl Il988 : 
Hubeny & Lan3 1 19921 1 19951 iLanz & Hubenyl l2003h for the 
model atmosphere calculations. We calculate plane-parallel 
model atmospheres in LTE taking into account bound-bound 
and bound-free transitions of ions listed in Table |2] We in- 
cluded lighter elements only. Although iron and other iron-peak 
elements belong to very important sources of opacity in the at- 
mospheres of hot s tars, we neglecte d their influence because 
iron is, according to iKhokhlova et al.l (12000.) . underabundant in 
the atmosphere of HD 37776 (see also Sect.©. Because it turns 
out that the bound-free transitions are important for the vari- 
ability of HD 37776, we note that the photoionization cross- 
sections of lighter elements in the TLUSTY ionic models are 
taken from the Opacity Project ( Seaton et alJl992 VButler et alJ 
19931 and references therein). Because here we study mainly 
the effect of hehum and silicon, for these elements we used 
the same models as used in the grid of model atmospheres ap- 
propriate to B-type stars calculated by .Lanz & Hubeny ( 20061) . 
Simpler atomic models were used for other elements to speed 
up the calculation of model atmospheres. 

For the spectrum synthesis (from which we calculate 
the photometric colours) we used SYNSPEC code. We took 
into account the bound-bound and bound-free transitions for 
the same ions as for the model atmosphere calculation (see 
TableEJ. 

We divided the visible star's surface into 90 x 360 ele- 
ments in latitude and longitude, respectively. We selected such 
a fine grid just from the numerical reasons to obtain a fairly 
smooth light curve. The local abundance of helium and sili- 
con (i.e., their number density relatively to the hydrogen) in 
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Fig. 1. The variation o f the hehum (left panel) and silicon abundances on the visible disc for different rotational phases (after 
Khokhlovaet"aDl2000l) 



each surface element w as set in accordance with the maps of 
Khokhlova et al. I (12000'). The emergent flux from these ele- 
ments was derived by the interpolation from the model grid. 
This grid of the model atmospheres contains 80 individual 
models. The grid is constructed in such a way that for each of 
the 10 He abundances, 8 different Si abundances were taken 
according to Table [3] The abundance of other elements for 
simp licity is assumed to be the solar one ( IGrevesse & Sauval 
Il998i) . i.e., the abundance of elements other than hehum and 
silicon relatively to hydrogen was kept fixed. Note that the ef- 
fective temperature and surface gravity of all these model at- 
mospheres are set to be the same. 

We calculate Eddington fluxes H(A, Y,Z) for model atmo- 
spheres with helium and silicon abundances Y, Z from the grid 
(A is the wavelength). The flux Hc{Y, Z) in a colour c is given by 
the convolution of H(A, Y, Z) with the transmissivity function of 
a given filter c. The transmissivity function is approximated by 
a Gauss function peaked at the central wavelength of the colour 
Al- with dispersion cr^, (see Table|4]i. This is performed for each 
colour c of the uvby system. The radiative flux in a colour c 
from individual surface elements HdQ.) is obtained by means 
of interpolation to a proper value of the He and Si abundances 
of the surface element concerned. Consequently, the flux HdQ.) 
is a function of spherical coordinates Q on the stellar surface. 
The total radiative flux observed at the distance D from the star 
with radius Rt is calculated as the integral over all visible sur- 



Table 2. Ions and the number of their excitation states consid- 
ered for the calculation of HD 37776 model atmospheres 



Ion 


Levels 


Ion 


Levels 


Ion 


Levels 


Ion 


Levels 


Hi 


9 


Ni 


21 


Nei 


15 


Sin 


40 


Hii 


1 


Nil 


26 


Neil 


15 


Sim 


30 


He I 


24 


Nm 


32 


Nem 


14 


Si IV 


23 


Hen 


20 


Niv 


1 


Neiv 


1 


Siv 


1 


Hem 


1 


Oi 


12 


Mgi 


13 


Sii 


14 


Ci 


26 


On 


13 


Mg 11 


14 


Siii 


10 


Cii 


14 


Oiii 


29 


Mgiii 


14 


S IV 


15 


Cm 


12 


Oiv 


1 


Mg IV 


1 


Sv 


1 


Civ 


1 















Table 3. Hehum and silicon abundances of the model grid 

[He/H] ^6 ^5 ^ ^3 1 ~l 1 2 T 
[Si/H] -3-2-1012 2.5 3 



face elements (lMihalaslll978h 



^(1)7 



lc{0, Q) cosfldQ 



visible 
surface 



visible 
surface 



(Q)Mc(6')cos6»dn, (1) 



4 



J. Krticka et al.: The light variability of the helium strong star HD 37776 



Table 4. Central wavelengths and dispersions of the Gauss fil- 
ter simulating the transmissivity functions of uvb y (Stromgre n) 
photometric system. These values are taken from lCox (l2000h . 



colour 


u 


V 


b 


y 


^ [A] 


3500 


4100 


4700 


5500 


a-a [A] 


230 


120 


120 


120 



Table 5. Calculated limb darkening coeflicients (Eq. (O) in in- 
dividual colours 



colour 


u 


V 


b 


y 




0.233 


0.194 


0.184 


0.164 


be 


0.163 


0.216 


0.204 


0.177 



where 6 is the angle between normal to the surface element 
and line of sight, udQ) - Ic{0,^)/IciO,Q.) describes the limb 
darkening (assuming the same Umb darkening for each surface 
element), Ic is the intensity in colour c and 

f/, = — f M,(6») cos edQ. (2) 

47T J 

For our calculations we use a quadratic limb-darkening law 

Uc{6) = 1 - flc (1 - cos 0) - (1 - COS ef . (3) 

We adopted the same limb darkening regardless the actual 
abundance of the surface elements. Namely, the quadratic limb 
darkening coefficients Oc, be in each colour c were derived from 
the solar abundance model. Different methods are available for 



the de rivation of limb darkening coefficients (e.g. iGrv gar et al 



971. 



We selected a method similar to Wade & Rucinski 



(1 19851) . which conserves the total flux. This leads to a condi- 
tion 

where H(A) is the emergent Eddington flux and I(A, 6 = 0) 
is the emergent intensity at the 6 = ray. The value of the 
coefficient a is derived by the least-square fitting of Eq. (O 
with the condition Eq. (|4| to values of m,(/1) calculated using 
SYNSPEC for different angles For the least-square fitting 
we used weights equal to angles fij. This leads to a formula 



a = 



i://, (l -/^0(l - 2//,) [(6 - ^)(l -^if - 1 + Uj] 



(5) 



where we, for abbreviation, dropped the wavelength depen- 
dence of all quantities. The limb darkening coefficients in each 
colour flc, be were derived from Eqs. (HJi, (|5]) using the gaus- 
sian smoothing over the colour c with parameters taken from 
Table |4] Bound-bound transitions in the calculation of Uc were 
also taken in to account. The calculated linear limb darkening 
coefficients are given in Table |5] 

Finally, the observed magnitude difference is 



Anic = -2.5 log 



A 



(6) 



where fc is calculated from Eq. ([T) and /J*^^ is the reference 
flux obtained from the condition that the mean magnitude over 
the phase is zero. As the spectral energy distribution within the 
width of a Stromgren filter is affected by the int erstellar redden- 
ing on ly marginally, provided that e. g. after IShore & Brown 
(Il990t) E(B -V) = 0.081, we can neglect the influence of the 
interstellar reddening. Note that the interstellar reddening in- 
fluences mainly the difference between two colours (the colour 
index), not the value of the change in a given colour. 

The calculation of magnitudes after Eq. (|6]l is performed 
over 36 rotational phases. 

3. Influence of chemical peculiarity on emergent 
fluxes 

As a result of higher opacity in the models with enhanced abun- 
dance of either helium or silicon, the temperature in the outer 
parts of the atmosphere increases due to the backwarming. This 
can be seen in the plot of the dependence of the temperature on 
the Rosseland optical depth Tmss for atmospheres with different 
chemical compositions in Fig. |2] To identify the main source 
of the opacity responsible for the backwarming we also calcu- 
lated atmosphere models neglecting the contribution of He-Si 
lines. Apparently, the line transitions influence the atmosphere 
temperature only for Tjoss S 10"^, where the stellar atmosphere 
is basically transparent in continuum. Consequently, the main 
influence of the enhanced abundance of helium and silicon on 
the atmosphere temperature stratification (in the region where 
the photometric flux forms, i.e., for 0.1 < Tioss S 1) is due 
to bound-free transitions. Indeed, as demonstrated in Fig. [3] 
in the case of silicon overabundance, the bound-free transi- 
tions from the higher excited levels of Si n and Si iii signifi- 
cantly contribute to the continuum opacity in the wavelength 
region where most of the flux is radiated from the atmosphere 
of HD 37776, i.e., in the far UV region of the Balmer contin- 
uum. 

The radiative flux emerging from the model atmospheres 
depends on the silicon and helium abundances, as can be seen 
in Figs, in |5] For solar chemical composition, the Balmer jump 
is visible in the plot. With increasing silicon abundance, the 
opacity in the ultraviolet (UV) region below roughly 1600 A 
increases mainly due to bound-free silicon transitions (and 
partly also due to the bound-bound transitions, see Fig.|3]l. The 
flux from this UV region is redistributed towards longer wave- 
lengths, partly into the optical region. Consequently, the silicon 
overabundant regions are brighter than the silicon poor ones in 
the uvby bands. A very similar situation occurs also in the case 
of helium overabundance. However, since for [He/H]=2 the he- 
lium dominates in the stellar atmosphere, jumps due to ioniza- 
tion of neutral helium occur. Helium overabundant regions are 
also brighter in the uvby bands than those with lower helium 
abundance. Decreased helium and silicon abundances below 
the solar values outside the bright spots do not significantly 
influence the flux in the uvby colours. As a result, we expect 
that only bright spots occur on the stellar surface, not the dark 
ones that may seemingly originate only due to a contrast effect. 

There is a relatively good agreement between the shapes 
of the calculated and predicted spectral energy distribu- 
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Fig. 2. The dependence of the temperature on the Rosseland optical depth Tmss in the atmospheres with various chemical com- 
positions. Crosses denote the case of atmosphere with enhanced abundance of either siUcon or helium, however neglecting the 
opacity due to the line transitions of these elements. 
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Fig. 3. Different sources of the opacity at the point where Tross ~ 1 for the solar chemical composition (left panel) and for 
[Si/H]=2. In the silicon overabundant atmosphere the bound-free transitions due to higher excited levels of Si ii and Si iii become 
important sources of the opacity in the far UV region of the Balmer continuum. 



tion (for solar abundanc e) in the visible region (taken from 
Adelman & Pypeiill985l note that from the observations only 
fluxes expressed relatively to the flux at 5000 A are available). 
The observed fluxes are slightly higher than the predicted ones 
in the Lyman continuum just below the Balmer jump. This dif- 
ference might reflect the influence of the chemical peculiarity. 

To identify the main sources of the light variability, we 
compared the fluxes calculated using models with neglected 
line transitions with that ones calculated with inclusion of the 
line transitions. This comparison showed that the line transi- 
tions contribute to the light amplitude only by up to percents. 

The change of the fluxes manifests itself through the change 
of the apparent magnitude. Fig. |6l The change is nearly the 
same in the b and y colours, which is caused by the fact that 
these colours lie on the Wien part of the energy distribution 
function and no important ionization edges are located there. 
The different change in the u colour is connected with the be- 
haviour of the Balmer jump and the occurrence of the jump due 



to helium ionization. The change in the v colour may be also 
slightly influenced by the presence of the Balmer jump. 

The change of the chemical composition of the atmosphere 
also influences the limb darkening as shown in Fig.|2] However, 
near the disk center for cos > 0.6 the variations of the limb 
darkening with abundance are not so significant. Consequently, 
we did not include these changes into our modelling, because 
the approximation of the limb darkening does not have a sig- 
nificant influence on the calculated light curves. For example, 
out test showed that there is only a relatively small difference, 
up to a: 0.001 mag, between the light curves calculated using 
quadratic and linear limb darkening. 



4. Simulation of the light curve 

4.1. Helium 

Since helium is the dominant element in the majority of regions 
on the HD 37776 surface dKhokhlova et al.ii2000i) . we started 
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Fig. 4. The emergent flux from the atmospheres with different helium abundances. The flux was smoothed by a Gaussian filter 
with a dispersion of 10 A to show the changes in continuum, which are important for photometric variability. The passbands of the 
uvby photometric system are also shown on the graph (gray a rea). The photometric reg ion was also enlarged and overplotted. For 
a reference we show also the observed mean flux derived by Adelman & Pyper ( 19851 , which we normalize using the predicted 
solar-abundance flux at 5000 A. 
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Fig. 5. The same as Fig.lH for silicon. 



studying the light variations caused by this element only. The 
corresponding light curves calculated taking into account the 
uneven distribution of helium on the HD 37776 surface (and 
flie solar chemical composition of other elements) are given 
on Fig. [8] Because the most helium-rich regions appear on the 
visible disk during the light minimum of the observed light 
variations, and because the helium overabundance induces the 
brightening of the stellar surface (Fig. |4|, the predicted light 
variations caused by helium are practically in antiphase with 
the observed ones. Consequently, although helium significantly 
contributes to the light variability of HD 37776, there has to 
also be another mechanism that dominates the optical variabil- 
ity of this star. 

To check the response of the light curve to the various val- 
ues of the abundance of helium we reduced it by 1 dex relative 
to the value derived by .Khokhlova et al.. (.2000,) . This, however. 



did not induce a remarkable change, as seen in Fig. [8] This is 
due to the fact that even with [He/H] = l the atmosphere starts 
to be dominated by helium already. 

4.2. Silicon 

The fact that silicon is overabundant on the observed disc 
aroun d the phase (j> = 0.75, when the star has its light maxi- 
mum Ad elman 1997, see also Fig. |9] and Fig.[TJ, along with 
our finding that silicon rich regions are bright, lead us to de- 
duce that silicon may be the main cause of the HD 37776 vari- 
ability. To test this, we calculated the light curve of HD 37776 
taking into account on ly the uneven distribution of silicon 
(iKhokhlova et al. 2000l) . The resulting light curve is displayed 
in Fig. |9] Apparently, taking into account only the uneven sur- 
face distribution of silicon, the amplitude of the light curve is 



J. Krticka et al.: The light variabiUty of the heUum strong star HD TillK) 



1 



-0.35 



-0.15 
-0.1 
-0.05 \ 


0.05 







He 


ich 


[He/H]=2 








solar metallicity 




/ \ 




He poor [He/H]=-2 


- 














u 


V 




b 




. Y. . . . 



CO 



3000 3500 4000 4500 5000 5500 6000 
^[A] 



-0.35 

-0.3 
-0.25 

-0.2 
-0.15 

-0.1 
-0.05 \ 


0.05 



Si rich [Si/H]=2 
solar metallicity 
Si poor [Si/H]=-2 



u 



3000 3500 4000 4500 5000 



5500 6000 



Fig. 6. The magnitude difference between the flux at the given abundance of helium {left panel) and silicon {right panel) H{A, Y, Z) 
and the flux for the solar abundance H{A, Yq,Zq) calculated as Am,i - -2.5 log{H{A, Y,Z)/H{A, Yq,Zq)). The passbands of the 
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too large when compared with the observed one (especially in 
the u colour). However, there is a good agreement between the 
light curve shapes. Note also that the amplitude is relatively 
sensitive to the value of [Si/H]. 

4.3. Helium and silicon 

The promising correspondence between the observed and pre- 
dicted light curves derived above encourages us to include both 
silicon and helium in the calculation of the light curve. To visu- 
alize the origin of the light variations, we plot the photometric 
spots on the surface of HD 37776 in different rotational phases 
in Fig.[TO]in the u colour. As it follows from this plot the shape 
of the light curve significantly depends on the inclination /. 

Inspecting Figs. |6l [8] we can conclude that the varia- 
tions due to helium and silicon are basically in the antiphase. 
Consequently, helium may suppress the too large amplitude de- 
rived assuming silicon surface variations only (especially in the 
u colour see Fig. |9|l. Indeed, the computed light curve shown 
in Fig.[TT] where both silicon and helium abundance variations 
are taken into account, confirms these considerations. There is a 
very good agreement between the predicted and observed light 
curves in the shape and amplitude. There is just a relatively 
small disagreement in the u colour. Some part of the remaining 
discrepancy between the predicted and observed curves may be 
due to uncertainties of spots location as derived from observa- 
tions, as the observed maps are tabulated at relatively coarse 
interval of 18° in longitude. 

Consequently, we conclude that the observed light variabil- 
ity of HD 37776 can be successfully simulated assuming the 
uneven surface distribution of individual elements, mostly sili- 
con and helium. 

5. Discussion 

The role of other elements An uneven surface distribution 
of elements other than silicon and helium may influence the 



light curve. Whereas the influence of oxyge n and iron, which 



are underabundant on the HD 37776 surface dKhokhlova et al 



is likely to be negligible, other elements (especially 
carbon for which the spectral line variability is observed, see 
Shore & Brown 1 19901) may modify the predicted light curve. 
On the other hand, if the iron abundance is much higher 
than what we assumed, then the iron opacity could be impor- 
tant for the spectral energy distribution and variability (pro- 
vided the uneve n surfa ce distribution of iron). We note that 
Khokhlova et al. I (l2000t) derived from Doppler-Zeeman map- 
ping that iron is relatively most abundant in spots where its 
deficit reaches up to 2 dex, while in the rest of the atmosphere 
the deficit is 4 - 5 dex relative to the solar value. These val- 
ues are surprisingly low and deserve further investigation using 
model atmospheres with magnetic fields included. 

NLTE effects The NLTE effects have a significant influence 
on the spectra of hot stars. These effects influence mainly the 
state of the outermost parts of the stellar atmospheres, where 



iiry 



ble continuum. Since lKhokhlova et al.l (120001) used LTE models 
for mapping of the HD 37776 surface (and, consequently, the 
real abundances may be different), we use also the LTE mod- 
els here. The neglect of the NLTE effects for both the abun- 
dance determination and synthetic spectrum calculation may 
be one of the reasons of the remaining small discrepancy in the 
u colour. 

Influence of the magnetic field Staude (1972) and Trasco 
(1972) showed that a strong surface magnetic field may in- 
fluence the emergent spectral energy distribution which may 
in photometry amount to hundredths of magnitude. Khan & 
Shulyak (2006) and Kochukhov et al. (2005) studied the influ- 
ence of the polarized radiative transfer and magnetic line blan- 
keting on the energy distribution, atmospheric structure and 
photometric colours and showed the magnetic field has a clear 
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Fig. 7. The continuum limb darkening at the wavelength corresponding to the center of the individual filter (Tab. HI in the model 
atmospheres with different chemical composition. The solid black line denotes the quadratic limb darkening (see Eq. (O, taking 
into account also lines) adopted in this paper 



relation to the visual flux depression. IStepienId 1978 b deduced 
that the eff'ect of the magnetic field on the stellar atmosphere 
leads to variation of T^s and log g over the stellar surface re- 
sulting in light variations. Consequently, the surface variations 
of the magnetic field may also contribute to the light variability. 

Vertical distribution of elements The vertically dependent 
abundance variations in chemical ly peculiar stars are suggested 
by many authors (e.g.. lRyabchikova.2004) . This effect may 
also influence our results, because the effective depth at which 
the lines form (i.e., at which the abundance map is in fact 
derived) differs from that at which continuum forms (i.e., at 
which the light variations are formed). Moreover, the chemical 
stratification itself may influence the emergent continuum flux. 
However, to our knowledge, no evidence of the vertical chem- 
ical stratification is available for HD 37776 atmosphere at the 
present time. 

Uncertainties of stellar parameters Uncertainties of the de- 
termination of the stellar parameters may influence the com- 
parison between observed and predicted quantities. However, 
since we study mainly differential flux variations in individ- 



ual colours, the influence of the determination of, e.g., the sur- 
face gravity acceleration or the stellar effective temperature 
are likely to be less important. Indeed, the comparison of the 
the magnitude differences calculated for [Si/H] -2 after Eq. |6] 
for stellar model calculated with the effective temperature by 
lOOOK lower or the surface gravity by 0.5 dex higher shows 
that the magnitude differences vary by about 0.02 mag (com- 
pare with the change of 0.1 - 0.35 mag showed in Fig.|6l). The 
magnitude differences are even less sensitive to the changes 
of the turbulent velocity, because the increase of the turbulent 
velocity by 2kms"' causes the variation of the magnitude dif- 
ference by only about 0.001 mag. 



Alternative surface chemical composition maps Because 
with the helium-only overabundance model we came up with 
a big discrepancy between the computed and observed light 
curves (Fig. [8]), we used another available model derived for 



the sur face chemical composition of helium only, bv lVeto et al 



(119911) . However, neither with this different surface pattern we 
were able to reproduce the observed light curve (cf. Fig. [8]l. 
A further map of the surface chemical composition which can 
be used for the rehable simulation of the fight curve is that 
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Fig. 8. The simulated light variations of HD 37776 in dif- 
ferent colours derived taking into account only the uneven 
distribution of helium on its su rface (solid Une) cornpared 
with individual uvby observations ( Pedersen & Thomsenll977 ; 
Adelman &Pvpedll985l;lAdelmanlll997l) . The dotted line de- 
notes the simulated light curve derived taking into account the 
helium abundance reduced by 1 dex. All observed magnitudes 
are shifted in such a way that their mean value is zero. A ver- 
tical offset of 0.05 mag was applied between each consecutive 
filter. 



of Bohlender & Landstreetl ( 1988 ). To test the influence of in- 
dependent surface chemical composition ma p on our result, 
we used the original map derived by Khokhlo va et al. I (SOOO), 
however the maximum and minimu m abunda nces of helium 
and silicon were taken from Bohlender & Landstreetl ([l988j). 
Generally, as a result o f lower silicon abunda nce in the spot 
(compared with that bv iKhokhlova et al. the amplitude 

due to silicon is lower when using .Bohlender & Landstreet 
(Il988h silicon abundance. However, this is compensated by 
lower maximum helium abundance. Consequently, our tests 
showed that also with this diff'erent surface map we are able 
to obtain a fair agreement with observations. 
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Adelman & Pyper (1985) 
Adelman (1997) 
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Fig. 9. The same as Fig. [8] taking into account only the un- 
even distribution of silicon. The dotted line denotes the light 
curve calculated for the surface distribution of silicon derived 
from the observations bv Khokhlova et al. (2000) decreased by 
0.5 dex. A vertical off'set of 0.05 mag was applied between each 
consecutive filter 



Comparison with observed flux distribution The next step in 
the comparison of the observed and predicted characteristics 
w ould be to compare the me an observed average flux derived 
by lAdelman & fVpen (Il985h and the theoretical one. Whereas 
there is a good agreement between these fluxes in the visible 
region, there is a slight disagreement (less than 10%) between 
them in the region between helium and hydrogen ionization 
edges in the u colour. The enhancement of the mean predicted 
flux in this region is caused by the enhanced helium abundance. 
Consequently, this might indicate that the maximum helium 
ab undance on HD 37776 surface is lower than that reported 
bv IKhokhlova et al. I (l2000l) . We also note that, for example, the 
effective temperature of HD 37776 was obtained using atmo- 
sphere models with solar chemical composition, which are un- 
realistic in the case of helium dominated atmosphere. However, 
the detailed study of this problem is beyond the scope of the 
present paper, which is concerned with the flux variations in 
individual passbands and not in the differences between mag- 
nitudes of individual passbands. 



Light variations in UV If the proposed mechanism of the fight 
variability is correct, the studied star should be variable also in 
the UV region and the observed light curve in the far UV region 
should be in antiphase with the optical and near-UV light curve. 
Indeed, a similar behaviour was found for other CP stars CU 
Vir and 56 Aii by ,Sokolov, (.2000. .2006) . 



Light absorption (ext i nctiori ) in the circumstellar environ- 
ment iTownsend et al. (2005) were able to simulate success- 
fully the light curve of anot her helium-rich chemically pecu- 
liar s tar cr Ori E (see also iNakaiimal 1 19851 : ISmith & Groote 
'2001) usi ng the model of rigidly rotating magnetosphere 
(iTownsend & Owockil l2005h . However, we argue that our 
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Fig. 10. The location of the calculated photometric spots in the 
u colour. We plot the smoothed flux emergent from the indi- 
vidual surface elements calculated using the observed surface 
distrib ution of helium and silicon derived by Khokhlova et al. 
(I2OO0I) in magnitudes. Here we do not include the limb darken- 
ing. 



-0.2 0.2 0.4 0.6 0.8 1 1.2 
phase 

Fig. 11. The same as Fig. [8] taking into account the uneven sur- 
face distribution of both the heUum and siUcon (solid line). A 
vertical off'set of 0.05 mag was applied between each consecu- 
tive filter. 

that the anomalous Zeeman effect is capable to induce such an 
enhancement of the Si 11 4131 A without the need of the large 
overabundance of the element. 



model of a rotating spotted star is more convenient for 
HD 37776. First, our method allows the simulation of the light 
curve just from the surface chemical distribution and we do not 
use any free parameter to fit the amplitude. Second, the light 
curve of HD 37776 is quite different from that of cr Ori E, 
which shows deep and r elatively narrow minima w ith maxi- 
mum amplitude in the u dPedersen & Thomsen ll977). Finally, 
cr Ori E is known for its emission in Ha line ([Walborn 1974), 
which is absent in the case of HD 37776. The uneven dis- 
tribution of chemic al elements on the surface of cr Ori E 
( Reiners et al. 2000h may, however, also influence the light 



curve of cr Ori E and may explain some secondary effects 
which are not explained by the cloud absorption model. 

Light variations of otiier cliemically peculiar stars Since siU- 
con is overabundant in the majority of chemically peculiar stars 
that show light variations (Si, He -weak, and He- strong stars, 
very frequently cool CP stars, e.g., lAdelmanlll973h it may play 
a crucial role also in their variability. However, this has to be 
tested also for other CP stars, as, e. g. Carpenter (1985) showed 



6. Conclusions 

In this work we were able to explain consistently the light curve 
of He-strong chemically peculiar star HD 37776. The observed 
light variations are explained in terms of the uneven surface 
distribution of helium and silicon. The helium and silicon spots 
modify the emergent flux predominantly due to the bound-free 
transitions. The predicted light curves reproduce the observed 
ones very well in their overall shape and amplitude. Although 
the agreement between the observed and predicted light curves 
may seem satisfactory, we stress that there are also other effects 
(e.g., NLTE effects, influence of other chemical elements) that 
may be important for the light curve. 

We have shown that bound-free transitions of mainly heav- 
ier elements (in our case silicon) accompanied with uneven sur- 
face distribution of these elements are crucial for the light vari- 
ability of magnetic chemically peculiar stars. To our knowl- 
edge, this mechanism has never been invoked in the literature 
to predict the light variability. 

Our results are important not only from the point of view of 
explanation of light curves of chemically peculiar stars, but the 
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comparison of predicted and observed light curves represents 
an important check of model atmospheres and their predictions. 
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